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, Abstract 

We study the small dispersion limit for the Korteweg-de Vries (KdV) equation 
u t + 6uu x + e 2 u xxx = in a critical scaling regime where x approaches the trail- 
ing edge of the region where the KdV solution shows oscillatory behavior. Using 
C*") ' the Riemann-Hilbert approach, we obtain an asymptotic expansion for the KdV 

solution in a double scaling limit, which shows that the oscillations degenerate to 
c j ' sharp pulses near the trailing edge. Locally those pulses resemble soliton solutions 

Qh. of the KdV equation. 

J: 
^__) . 

1 Introduction 

We consider the Cauchy problem for the Korteweg-de Vries (KdV) equation 

ut + 6uu x + e 2 u xxx = 0, u(x, t = 0, e) = uq(x), e > 0, t > 0, (1-1) 

: 

in the small dispersion limit where e — > 0. We consider real analytic negative initial 
data with sufficient decay at infinity and with a single negative hump. For small t > 0, 
the solution to this problem can be approximated [211 E] by the solution to the Cauchy 
problem for the (dispersionless) Hopf equation 
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ut + 6uu x = 0, u(x, 0) = uo(x), t > 0, (1-2) 

which can be solved using the method of characteristics. At time 
1 



max ?6K [-6'a' (^)]' 

the Hopf equation reaches a point of gradient catastrophe where the derivative of the 
Hopf solution blows up. For t > t c , the Hopf solution only has a multi- valued con- 
tinuation. However, the KdV solution is well-defined for all positive t and e > 0: the 
dispersive term e 2 u xxx regularizes the gradient catastrophe. For t slightly bigger than 
t c , the KdV solution u{x, t, e) develops rapid oscillations |21|, 124] that in the limit e — > 
are confined in a certain interval [x~ (t), x + (t)}, see Figure [TJ In the (x,t)-plane, the 
oscillations take place in a cusp-shaped region (which depends on the initial data) as 
illustrated in Figure 2. Inside the cusp-shaped region for t slightly bigger than t c , the 
KdV solution u(x, t, e) can be written, asymptotically for small e, in terms of the Jacobi 
elliptic theta function and the complete elliptic integrals of the first and second kind 
E(s) and K(s) [2H EJ Ml CE] : 

d 2 

u(x,t,e) ~^i + /3 2 + ^ 3 + 2a + 2e 2 ^lntf(ft(x,t);T), (1.3) 
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Figure 1: Solution of the KdV equation with e = 10 2 , initial data Uq(x) = — sech 2 (x), 
and t = 0.4. 

where Q, a, and T have the form 



n(x, t) = V 2 ^ ( ^ 3 [x - 2t{0 x + 2 + 3 ) - q], (1.4) 
a(s) = + (A -ft)-^7-r, T = *WT> s = « a"' C 1 - 5 ) 



Note that K'(s) = K{yX — s 2 ), and $ is defined by the Fourier series 
tf(z;T)=]Te* 



„7nn 2 T +2winz 
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The formula for q in the phase Q in (11. 4h is equal to pU [8] 

A ' A) = 27^ 7-! 7-! ^ v^Vl-^ 2 ' (L6) 

where /z,(y) is the inverse function of the decreasing part of the initial data uq. Formula 
(|1.3p can be written also in terms of Jacobi elliptic function dn in the form 

u(x,t,e) ~/? 2 + /3 3 -/?i + 2 - f 1_/?2 



dn\2K(s)Q) (1.7) 
02 + fo-Pi + 2{0i - /3 3 )dn 2 (2K(s)0 + (2k + l)K(s)), 

for any integer k since dn(it; s) is periodic with period 2K(s). 
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Figure 2: Cusp-shaped oscillatory region in the (x, f)-plane for t > t c . 



For constant values of the /3j's, the right hand side of (II. 7|) is an exact solution of 
KdV. However in the description of the leading order asymptotics of u(x, t, e) as e — > 0, 
the numbers (3% > fa > /?3 depend on x and t and evolve according to the Whitham 
equations [27] 
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with a as in (|1.5p . 

The Whitham equations (jl .81) can be integrated through the so-called hodograph 
transform, which generalizes the method of characteristics, and which gives the solution 
in the implicit form [25] 



1,2,3, 



(1.9) 



where the v^s are defined in (|1.8p and = Wi{(3\, fa, /%) for i = 1,2,3 is obtained 
from an algebro-geometric procedure by the formula 




fc=i 



dq 

Wi 



+ q, 



1,2,3, 



(1.10) 



with q defined in (jl.6p . The formula (jl.6p for ^ is valid as long as (3% does not reach 
the minimal value of the initial data no. When reaches the negative hump it is 
also necessary to take into account the increasing part of the initial data fn- However, 
formula (11.61) is sufficient for the purpose of this manuscript. 



Near the boundary of the oscillatory cusp-shaped region, neither the Hopf asymp- 
totics nor the elliptic asymptotics are satisfactory. Three different transitional regimes 
can be distinguished: (1) the cusp point where the gradient catastrophe for the Hopf 
equation takes place and where Pi = = @3, (2) the leading edge of the oscillatory 
zone where = @3 , and (3) the trailing edge of the oscillatory zone where /?i = /?2 • 

Near the point of gradient catastrophe, it was conjectured pUlE] and proved after- 
wards [5] that the asymptotics for the KdV solution are given in terms of a distinguished 
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Painleve transcendent, namely a special smooth solution U(X,T) to the fourth order 
ODE 



X = TU 



qU 3 + -(U 2 X + 2UU XX ) + m U xxxx 



which is the second member of the first Painleve hierarchy. In a double scaling limit 
where e — ► and simultaneously x and t approach the point and time of gradient 
catastrophe x c and t c at an appropriate rate, the KdV solution has an expansion of the 
following form [5], 

«(*, t,e)=u c +( 2 4) V7 U ( X " Xc " 6U f " tc) , ^4 ) + O (e^ 



k 2 



(8fce 6 )7 (4fc 3 e 4 



where A: is a constant depending on the initial data. This expansion holds for negative 
real analytic initial data with a single negative hump and sufficient decay at infinity, 
but it is conjectured by Dubrovin that it has a universal nature and extends to all 
Hamiltonian perturbations of hyperbolic systems near the point of gradient catastrophe. 
It is remarkable that the function U (X, T) itself is a solution to the (re-scaled) KdV 
equation Ut + UU X + ^Uxxx = 0. 

Near the leading edge at the left of the zone where the oscillations appear (see Figure 
[T]) , numerical results [T7] showed that the amplitude of the oscillations is asymptotically 
described by the Hasting-McLeod solution to the second Painleve equation. For t 
slightly bigger than t c , in a double scaling limit where e — * and simultaneously x 
approaches the leading edge x~ at an appropriate rate, it was proved in [6] that 



4e l/3 

u(x, t,e) = u T-prA 



x — x 



cos (^)+o( e n 



c 1 / 3 \Ju — v e 2 / 3 _ 

with A the Hasting-McLeod solution to the second Painleve equation, the phase 
G(x,t) = 2v^^(x-x-) + 2 /"(/!(£) + 6*) V? 3 ^, 

J V 

and x~(t), u(t) and v(t) are obtained from (|1.9p in the limit /?2 = @3 = v. This 
expansion holds, just like near the gradient catastrophe, for negative real analytic 
initial data with a single negative hump and sufficient decay at infinity. It is natural 
to expect that it is, to a certain extent, universal for a whole class of equations. 

Our aim is to give an asymptotic description of the KdV solution near the trailing 
edge. Here the behavior of the KdV solutions is genuinely different compared to the 
Painleve type behaviors near the point of gradient catastrophe and the leading edge. 
This is not surprising when considering Figure [H where we observe that the amplitude 
of the oscillations is of order 0(1) near the trailing edge, whereas it smoothly decays 
towards the leading edge. We will show that, in the limit where e — > 0, the last 
oscillations at the right end of the oscillatory zone behave like solitons which are, in 
the local scale, at a large distance away from each other. Recall that the KdV equation 
admits soliton solutions of the form a sech 2 (bx — ct) . We should note that the trailing 
edge asymptotics we will obtain show remarkable similarities with recently obtained 
critical asymptotics for the focusing nonlinear Schrddinger equation [3]. 
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The trailing edge x + {t) of the oscillatory interval (i.e. the right edge of the cusp- 
shaped region in Figure 2) is uniquely determined by the equations [Ml IE] 



x + (t) = 6tu(t) + f L (u{t)), 
6t + 6(v(t);u{t)) = 0, 
v(t) 



lu(t) 

Here v (t) > u(t) and 



(6* + 0(A; u(t)))y/X-u(t)dX = 0. 



(1.11) 
(1.12) 

(1.13) 



e(v;u) 



1 



2^1 



(1.14) 



and /l is the inverse of the decreasing part of uq. System (jl.lip - (|1.13p is the confluent 
form of the Whitham equations where /3i = /?2- 

Before stating our result, let us take a look at the elliptic asymptotics for KdV in 
the limit (3\ — > /?2- Note that this is a formal limit of the asymptotic expansion for 
KdV and that there is no rigorous argument to justify that this actually leads to a 
good approximation for the KdV solution near the trailing edge. The phase defined in 
(fl~4"j) behaves like 



x — x 1 , 

~ yw — u, 



as f}\ — fa — > 0, 



and 



K(s) ~ - In 
V ; 2 



1 - s 



as s — > 1. 



In addition the Jacobi elliptic function dn(u; s) — > sech(u) as the modulus s — ► 1. 
Therefore the formal limit of the elliptic solution (jl.7p as s — ► 1, /3i, 02 — * v, /?3 — ► u, 
gives 



n(x, t, e) ~ n + 2(u — u) sech 2 



1, 



■y/v — U + (k + -) In 



1 



(1.15) 



We would like to underline that the limit we computed in the above expression has 
appeared many times in the literature, starting from the seminal paper of Gurevich- 
Pitaevski [19] but the important second term of the phase in (|1.15p which was calculated 
later by Deift, Venakides and Zhou [8], was and is still ignored in some literature. When 
taking the limit s — > 1 for fixed e, the solitonic sech 2 -term in (|1.15p is of order 0(fi\— fyi)- 
However it is clear from Figure [T] that the amplitude of the oscillations near x + is of 
order 0(1). This indicates that it is necessary to perform a double scaling limit letting 
e — > and x — > x + simultaneously at an appropriate rate. If we want to capture the 
top of an oscillation close to the trailing edge, we should take the double scaling limit 
in such a way that the phase of the sech 2 term is zero for some integer k. If k > 0, 
this will turn out to be consistent with the rigorous asymptotic expansion for u we will 
obtain below in Theorem 11.21 up to a phase shift. 

Similarly as in [U [6], we consider initial data u${x) which satisfy the following 
conditions. 
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Assumptions 1.1 

(a) Uq(x) is real analytic and has an analytic continuation to the complex plane in a 
domain of the form 

S = {z G C : |Imz| < tan6>o|Re2i|} U {z £ C : |Im z\ < a}, 

for some < 0q < tt/2 and a > 0, 

(b) uq(x) decays as x — > oo in S such that 

u (x) = O (ps+j) > « > 0, (1-16) 

(c) for real x, uq(x) is negative and it has a single local minimum at a certain point 
xm, with 

u' (xm) = 0, u'q(x m ) > 0, 
for simplicity we assume that uq is normalized such that uq{xm) = — 1? 

(d) the point of gradient catastrophe for the Hopf equation is 'generic' in the sense 
that 

fl' K) ^ 0, (1.17) 

where fz, is the inverse of the decreasing part of the initial data uq. 

Because of condition (b), we will be able to apply direct and inverse scattering theory, 
and to setup a Riemann-Hilbert (RH) problem for the Cauchy problem of KdV. Con- 
ditions (a) and (c) will enable us to keep control over the reflection coefficient in the 
small dispersion limit; condition (d) will be needed only in the core of the RH analysis, 
but is nevertheless needed for the result stated below. 

Theorem 1.2 Let uq{x) be initial data for the Cauchy problem of the KdV equation 
satisfying Assumptions and let x + = x + (t), u = u(t), and v = v(t) solve the 
system \1.11\) - [1.13\) . There exists T > t c such that for t c < t < T, we have the 
following expansion for the KdV solution u(x, t, e) as e — > 0, 

(e In \ °° 

x++ € ^= y,t,e) =u + 2(v-u) V sech 2 (X fc ) + C(eln 2 e), (1.18) 
J fo 



where 



11 1 
x k = -{--y + k)hxe- ln(V2nh k ) ~{k + ^) In 7, 

2 i . , (L19) 

hk = 1 -= , 7 = 4(i> -u)4Tj-d v 6(v;u), 
7T* v k\ 

and 6 is given by |i. 1J$ . This expansion holds point-wise for y £M and uniformly for 
y bounded. 
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Figure 3: Cartoon picture of u yx + + ^M= y : t, e) as a function of s = —y for e = 10 J , 
or in other words a magnified version of Figure [1] near the trailing edge with the x-axis 
stretched with a factor C(l/(e| In e|)). One observes narrow pulses near the negative 
half integers. 

Remark 1.3 Observe that each term in the sum of (11. 18ft generates a pulse with 
amplitude 2(v — u) for y near a half positive integer (as shown in Figure 3), which 
can be seen as a soliton. The term sech 2 (Xfc) is centered at y ~ k + i, and it already 

decreased to order 0{e^) for y near k and k + 1. Near k— ~ and the contribution 

of sech 2 (A^,) is absorbed by the error term 0{e\v? e). For any y, the infinite sum in 
(j!.18j) thus reduces to the sum of the two solitons centered closest to y. The sum of all 
the other terms is of order 0(eln 2 e). Values of y for which y < — ^ lead us out of the 
oscillatory zone: the contribution of the sum of solitons to (|1.18|) is small, in this case 
we have u(x, t, e) = u(t) + 0(eln 2 e). 

Remark 1.4 Whereas the function u is clearly decreasing for x to the right of the 
trailing edge in Figure [IJ this is not visible in Figure 3 or in the asymptotic expansion 
(jl,18p . The reason is that this effect is of order elne in the local variable y. On the 
other hand, it is not clearly visible in Figure Q] that the oscillations degenerate to sharp 
pulses towards the trailing edge. This is a consequence of the fact that the length of 
the pulses scales with lne|), which decreases slowly as e — > 0. 

Remark 1.5 Analyzing the Whitham equations carefully in the limit (3% — ► /?2, it 
is possible to deduce an expansion of the form (|1.18[) from (|1.15p . but the formally 
obtained misses a constant term (depending on k) compared to (|1.19p . 

Remark 1.6 The time T > t c appearing in Theorem 11.21 should be sufficiently small 
such that Proposition l3.1l holds. Loosely speaking, this means that the G-function which 
we will construct below is not allowed to have any singular behavior for t c < t < T. 

Remark 1.7 The techniques we will use to prove Theorem 11.21 are very close to the 
ones used in [U HI [22], where orthogonal polynomials on the real line with respect to 
a critical exponential weight were studied. Those polynomials describe the birth of a 
cut in unitary random matrix ensembles, see also [12]. Although the focus in [2"1 14^ 122] 
was on the random matrix eigenvalues rather than on the orthogonal polynomials, we 
expect that the recurrence coefficients for the associated orthogonal polynomials admit 
an asymptotic expansion similar to (|1.18p in an appropriate double scaling limit. 
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In Section [21 we will set up the RH problem for the KdV equation, and we will 
recall asymptotic results about the reflection coefficient in the small dispersion limit. 
In Section [3] we will analyze asymptotically the RH problem using the Deift/Zhou 
steepest descent method. At the heart of the analysis lies the construction of a local 
parametrix built out of Hermite polynomials. The degree of the Hermite polynomials 
depends on the scaling variable y. The transitions where the degree of the Hermite 
polynomials increases takes place for y near positive half integers and are responsible 
for the presence of the spikes in the asymptotic behavior of the KdV solution u. 

2 RH problem for KdV and reflection coefficient 

It is well-known that solutions to the KdV equation can be expressed in terms of a 
RH problem. This fact relies on the direct and inverse scattering transform [H I13j . 
Consider the following RH problem. 

2.1 RH problem for M 

(a) M(X;x,t,e) : C\K -> C 2x2 is analytic. 

(b) M has continuous boundary values M + (A) and M_(A) when approaching A G 
K \ {0} from above and below, and 

/ 1 (\ \ 2ia(\;x,t)/e\ 

M + (A)=M4A)(_ r _ (A;£)e _ 2ia(A;x , t)/e ( r, r(A;e) | 2 } forA<0, 
M + (A) = M_(A)<x 1; ax = ( J J J , for A > 0, 

with a given by 

a(X; x, t) = 4t(-A) 3/2 + x(-A) 1/2 , (2.1) 

defined with the branches of (— X) e which are analytic in C \ [0, +oo) and positive 
for A < 0. Furthermore for fixed x, t, e, we impose M(A) to be bounded near 0. 

(c) M(A;x,t,e) = (/ + C(A- 1 )) I i as A ^ oo. 



viv— A — iy — Xi 



In particular we can write 



Mn (A; x,t,e) = l+ Ml .^(M> e ) + {l/X), as A ^ oo. 

In general it is true that, if the RH problem for M is solvable in a neighborhood of xq 
and to, then the function 

u(x,t,e) := -2ied x M 1A1 (x,t,e), (2.2) 

is, locally near xq and to, a solution to the KdV equation. Modifying the function 
r(A; e) in the jump matrix leads to different KdV solutions. If r(A;e) is the reflection 
coefficient corresponding to the Schrodinger equation 

d 2 

^dx^ 1 = U ° (X)/ ' 
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(with uq satisfying Assumptions 11.1] but also under much milder conditions) then u is 
the solution to the Cauchy problem (jl.ip for KdV with initial data uq. In this case the 
RH problem is solvable for all x G M and t > 0; its solution can be constructed using 
fundamental solutions to the Schrodinger equation. 

Our main task in what follows, is to approximate the RH solution M asymptotically 
for e — * and x close to the trailing edge. To this end, we need estimates for the 
reflection coefficient r(A; e) as e — ► 0. 

2.2 Asymptotics for the reflection coefficient 

Semiclassical asymptotics for the reflection coefficient as e —* were obtained in |23} [T5] 
for initial data uq that are such that Assumptions 11.11 hold: we have 

r{\;e)=ie-T^ x \l + 0{e)), for A G (—1,0), (2.3) 

r(A;e) = 0(e"f), for A G (-oo, -1), (2.4) 

where p is given by 

MA) 



p(X) = f L (X)V^X + yu (x)-X-V^X)dx. (2.5) 



Moreover, there is a sector containing R such that r(A; e) is analytic for A in this 
sector 1151]. Let us write 



k(A; e) = -tr(A; e)e^ (A) , for A G SI+, (2.6) 

where 

0+ := {A G C : -1 -5 < ReA < -6,0 < ImA < 5, |A + 1| > -}. (2.7) 

For sufficiently small <5 > 0, we have the following as e — > 0, see [23j and [6, Section 2.2] 

k(A; e) = 1 + 0(e), uniformly for A G 0+, (2.8) 

r(A;e) = 0(e - ?), for A < -1 - 5, c> 0. (2.9) 

In addition 

1 - |r(A;e)| 2 ~ e"^ T(A) , for-l + 5<A<0, (2.10) 

with 

/k(A) 



t(A) = / - n (x)(ix, (2.11) 

and fa is the inverse of the increasing part of the initial data uo- Those estimates for 
the reflection coefficient are essential for the RH analysis we will perform in the next 
section. 
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3 Asymptotic analysis of the RH problem 



We study the RH problem for M asymptotically in the small dispersion limit when 
x is close to the trailing edge. Our approach is based on the Deift/Zhou steepest 
descent method [9] which has been applied to the KdV RH problem in [8]. We follow 
roughly the same lines as in [6], where the RH problem was studied near the point 
of gradient catastrophe and near the leading edge. Compared to the situation near the 
leading edge, the main difference here is the construction of the local parametrix near 
the point v, which will be built out of Hermite polynomials. 

3.1 G-function and transformation MhT 

Let us define, for t c < t < T, G = G(A; x, t) by 

Ju Vn - u kn - a) 

where u = u(t) is the solution to (|l.ll|) - ([1.13p . and a and p given by (|2.5|) and (|2.ip . 
The square root (u — A) 1 / 2 is analytic for A G C \ [u, +oo) and positive for A < u. G 
has the asymptotic behavior 

G(A; x) = Gi(x, t) . * + 0(A _1 ), as A - oo, (3.2) 

with 

G 1 (x,t) = l f° P(V) -airn_^t) ^ 



Since p(rj) is independent of x, we have 



d x G 1 (x,t) = -- [ J^L drj= u (3.4) 



7T Ju V V ~ U 2 

G is analytic for AgC\[u, +oo) and it satisfies the properties 

G+(A) + G_(A)-2 / 9(A) + 2a(A) = 0, asAe(u,0), (3.5) 

G+(A) + G_(A) = 0, as A G (0, +oo). (3.6) 

In order to modify the jumps for M in a convenient way without losing its normalization 
at infinity, we define 

T(A) = M(A)e-^ A;:c K (3.7) 
Then, using the RH conditions for M, we obtain the following RH problem for T. 

RH problem for T 

(a) T is analytic in C \ R, 
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Figure 4: The jump contour Eg after opening of the lens 



(b) T+(A) = T_(A)«t(A), as A G R, with 



v T (\) 



1 



r(A)e" 



KG(A)+a(A)) s 



v _ f (A) e -T(GW+«(A)) i_| r (A)| 2 /' 

V;( G +W- G -W) r(A)eT^ A ) 

-f(X)e-TPW (l-|r(A)| 2 )ei(°+W- G -W) 



as A > 0, 

as A G (— oo, u), 

as A G («, 0), 



(3i 



(c) T(A) = (/ + 0(A^))(.^ 



-A 



as A —* oo. 



(3.9) 



(3.10) 



By (H3D, (|33D, and ([57T|> . we obtain the identity 

u(x, t,e) = u — 2ied x Ti ) u(x, t, e), 
where 

r n (A;x,t,e) = l + ^i^^ + 0(A- 1 ), as A ^ oo. 

3.2 Opening of lenses T \— > S 

Let us define the function <j) by 

4>(\; x, t) = G(A; x, t) — p(X) + a(A; x, i), 

such that (j> is analytic in a neighborhood of (— l,u), with an analytic extension to the 
whole region f2 + defined in (|2.7p . Under the condition that (jl.ll|) holds, it was shown 
in Lemma 3.2] that 



0(a ; x, t) = v^a(x - x + ) + (/i(e) + et) vf^a^, 

with, as usual, principal branches of the square roots. It then follows that 



(3.11) 



4>'{\;x,t) 



l 

2\/u - A 



(x — x ) - V« - A[6i + 0(A; u)], 



(3.12) 
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where 9(X;u) is given by (|1.14p . In a complex neighborhood of v it is convenient to 
define <j> as the analytic continuation of — i(f> from the upper half plane, 

0(A; x, t) = =F^(A; x, t) 

rX 

= -y/X - u(x -x + )+ / (f L {£) + 6t)y/\-£d£, as ±ImA > 0. (3.13) 

J u 

We then have using (fLl^ - (fl~T3l) . 

0(v;x + ,t) = 0, 4>'(v;x + ,t) = 0, 4>"(v; x + , t) = y/v^u 6'(v; u) < 0. (3.14) 
Now we can express the jump matrix for T in terms of (ft and 4>: 

v T {X)=\ 1 ' ' 22r,A, asXe(u,0), (3.15) 

V ' \iK*(X;e) (l-|r(A)| 2 )e-!^ A )y K , J, \ J 

( 1 m(A;e)e^ A A . r , . , ,„ 1<s . 



where (f>*(X) = 4>(X), and 0(A) and 4>*(X) should be understood as the + boundary 
values on (—1 — 5,-1). 

Similarly as in [6, Proposition 3.1], <j) and <f> satisfy a number of inequalities. 

Proposition 3.1 There exists a time T > t c such that for t c < t < T, we have 
4>{X;x + )<0, as X G (u,0\ \{v}, 

<f>'(X; x + ,t) > 0, as X G [-1, u), (3.17) 

-r(A) - (j)(X;x + ,t) < 0, as X G(u,0]. 

Proof. If (|1.17p holds, the function 6t c +9(X; u(t c )) has a double zero at A = u(t c ) and 
is strictly negative elsewhere on [—1,0], see [5J. For t slightly bigger than t c , because 
of the smoothness of u(t), 6t + 9(X;u(t)) can have at most two zeroes on [—1,0]. It 
follows from (|1.12|) that one zero lies at A = v(t). Because of ()1.13p . another zero 
must lie in between u and v. Consequently we have that 6t + 9(X;u(t)) is negative for 
A G [— l,u) U (v,0] for sufficiently small times after the time of gradient catastrophe. 
Now it is straightforward to verify the first and the second inequaltity. 

For the last inequality, it is straightforward to verify by (|2.1ip that 

-r(A) - 0(A;x") = -4t(A - u)§ + j" y/X^f' L (0^ ~ £ 7^/^(0^, 

where fa is the inverse function of the increasing part of the initial data Uq(x). The 
inequality follows directly because each of the three terms on the right hand side is 
negative. □ 

On the interval (—1 — 5, u), we can factorize vt- 
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Using this factorization, we can open lenses as shown in Figure 13.21 We choose the 
lenses such that the upper lens Si lies in Q + (for some fixed but sufficiently small 
6 > 0) and such that the lower lens E 2 is the complex conjugate of £i. Define 



5(A) 



T(A) 



T(A) 



1 -m(A;e)eT^) N 







1 J 

1 o N 

2i, 



in the upper part of the lens, 



>, , . , , . in the lower part of the lens, 

T(A), elsewhere. 
By the RH problem for T, we obtain modified RH conditions for 5. 

RH problem for 5 

(a) 5 is analytic in C \ £5, 

(b) 5+(A) = 5_(A)u s for A € E s , with 



,0 1 J' 

1 o\ 

v m*(A;e)e-T^(^) lj' 

'efft**) i*(A;e) 
,m*(A;e) (1- |r(A)| 2 )e"7^) 



on Si, 

on E 2 , 

, as A G (it, 0), 
as A G (0, +00), 



as A G (— 00, —1 — 5). 



(c) S(\) = (I + 0(\-i))(.J- 



1 



-A -W-X 

Since T(A) = 5(A) for large A, we have, using (|3,9 

u(x, t,e) = u — 2ied x S% ! ii(x, t, e), 
where 

5n(A;x,t,6) = l+ gl ' 1 ^_ t ' £) +0(A- 1 ), 



as A — > 00. 



as A — > 00. 



(3.19) 



(3.20) 



(3.21) 



The construction of 5 has been organized in such a way that the jump matrices for 
5 decay uniformly to constant matrices, except in arbitrary small neighborhoods of u 
and v. 

Proposition 3.2 There exists T > t c such that the following holds for t c < t < T. For 
any fixed neighborhoods U u ofu and U v ofv, there exists 5q > such that for \x — x + \ < 
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<5o, the jump matrices vs{X) = (I + 0(e))v^ 00 ' (X; x, e) uniformly on \ (U u UU V \J {0}) 
if e — ► 0, with 



W («0(A) = <{ 



(7i, on (0, +00), 

io"i, on (n, 0), (3.22) 
J, elsewhere. 



Proof. See Proposition 3.3]. □ 
For x > x + + „ e ! n<: , it is straightforward to verify that the convergence is also uniform 

2yu— u 10 * 

onln[/„ for sufficiently small U v . 
3.3 Outside parametrix 

Let us first ignore small neighborhoods of u and v, and uniformly small jumps. Then 
we need to solve the following RH problem for the outside parametrix. 

RH problem for p(°°) 

(a) P(°°) : C \ [u, +00) -> C 2x2 is analytic, 

(b) satisfies the jump conditions 

p M = pW^^ on ( 0)+oo)) (3.23) 

pjoo) = ^pM^ on ( U;0 ) \ |^| 5 (3 24) 



(c) p(°°) has the following behavior as A 



00, 



PM(A) = (/ + 0(A- 1 )) ( . ( _* )1/2 (3.25) 



This RH problem is solved by 

P (oo) (A) = (-A) 1/4 (n - \)~ a ^ 4 N, N=(^ 
This solution is bounded near v. Depending on the value of 



I 1 

i —i 



(3.26) 



y:=2Vv~=lL^—, (3.27) 
erne 

we will need an outside parametrix which has singular behavior near v. Therefore we 
observe that 

Pjf°\\) = (-A) x / 4 (n- X)-^/ A ND{\)^, for±ImA>0, (3.28) 
is a solution to the RH problem for p(°°) for any k G N, with 



_ ✓p-^ (M9) 

V A — it + V" ~~ u 
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Indeed this follows from the fact that D is analytic on C \ (— oo, it], with 

D + (X)D_(X) = 1, forA<u, 
D(oo) = 1. 

Near v, we have that P^ oo) (A) A ±feT3 is bounded for ±ImA > 0. 



(3.30) 
(3.31) 



The appropriate choice for k turns out to be as follows: if y < 0, we take k = 0, 
and if y > 0, we let k be the non-negative integer closest to y (for the half integers, we 
may choose k = y — \ or k = y + i), so y — | < k < y + tj. As A — ► oo, we have 



P fc (oo) (A) (J + ^ + cXA 



3.4 Local parametrix near it 



1 1 

i(-A)V2 -i(-A)V2 

2ik\/v — u 



(-A)V 



a 3 + 0(A- 1 ) . (3.32) 



In [6] Section 3.5], the Airy function was used to construct a local parametrix in a 
neighborhood f/ u of the point u. If conditions (jl.lip - (|1.12p hold, it was proved that 
the constructed parametrix satisfies the following properties. 



RH problem for P u 

(a) P u is analytic in U u \ Eg. 

(b) P u satisfies the jump conditions 

1 ie~ 



P U ,+(X) = P U ,_(A) 
P U ,+(A) = P U ,-(A) 
P„,+(A) = P U ,_(A) 



^0(A;x) 



v 1 
1 







L 0*(A;x) j I ' 



e ^i0+(A;x) • 



as A e Ei n U u , 



as A G E 2 n J7 U , 



as A > it. 



(3.33) 



(3.34) 



(3.35) 



(c) As e — > and x — > x + in such a way that x — x + = 0(e 2//3 ), P u matches with 
p(°°) in the following way, 



P U (X)=P^(X) 



is(X;x) 2 s(X;x) 

;&3 — — <y\ + U{e t 



4e(-C(A))V2 



4C 



Here s(A; x) is an analytic function of A G U u with 
s(«;x) = (-f' L (u)-6ty 1/3 (x-x+). 



for A G dU u - 
(3.36) 

(3.37) 



In the (stronger) double scaling limit where e — > and x — > x + in such a way that 
x — x + = O(elne), it follows that 



P U (A) = P(°°)(A) [I + 0(eln 2 e)] , for A G 017 U . 



(3.38) 
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For the explicit construction of P u in terms of the Airy function, we refer to [51 Section 
3.5]. It will be important in what follows that, in the limit e — > 0, the jump matrices 
for P u are equal to those for S up to an error of order 0(e). This follows from f|3. 19j) 
and the asymptotic formulas (|2.8p and (|2.10p for the reflection coefficient. 

3.5 Local parametrix near v 

The construction of the local parametrix will depend on the value of 

y = 2^AT^ X ~ X . (3.39) 
erne 

If y < —1, there is no need to construct a separate local parametrix near v since the 
jump matrix for S is equal to the one for p(°°) up to an 0(eln 2 e) error; in this case 
we write for notational convenience P v (\) = I. Let us assume now that y > — 1. The 
aim of this section is to construct a local parametrix in a small neighborhood U v of 
v having approximately the same jump property as S has near v (in the limit e — > 
0), and matching with the outside parametrix at dU v . Substituting the semiclassical 
asymptotics (j2.8j) and (|2.10p for the reflection coefficient, the jump matrix for S given 
in (I3.19|) behaves as follows when e — > 0, 

v s (X) = (I + 0(e))l ee . M . (3.40) 

This brings us to the RH problem for the local parametrix. 

RH problem for P v 

(a) P is analytic in U v \ M, 

(b) P has the jump condition 

P + (X) = P-(X) I 6 . IJ, (3.41) 

(c) if we let e — > and simultaneously x — > x + in such a way that 

Ix-x+l < Me|lne|, (3.42) 
we have 

P(A) = p(°°)(A)(l + o(l)), asAeS^. (3.43) 

We will construct P explicitly in terms of a model RH problem built out of Hermite 
polynomials. 

3.5.1 Model RH problem 

Define by 



-2mhk-\Hk-i{Q 



h k K ^> 2ivih k JR u-C I C 2 

*(C;fc)=| U-^ 3 , 




for C € C \ R, (3.44) 
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where denotes the degree k Hermite polynomial, orthonormal with respect to the 

2 

weight e~ x on R. The leading coefficient of the normalized polynomial is equal to 

** = *W <3 - 45) 

and we agree H_\ = h-% = 0. = solves the following RH problem for 

k E N U {0}, which is a particular case of the RH problem for orthogonal polynomials 
discovered in [T5] . 

RH problem for ^ 

(a) * : C \ R -» C 2x2 is analytic. 

(b) has continuous boundary values for £ £ R, related by the condition 

*+(C) = *-(C)(j })■ (3-46) 

(c) As C — > oo, we have 

(3.47) 

The further terms in the large-£ expansion can be calculated explicitly, but are unim- 
portant for us. 

3.5.2 Modified model RH problem 

In order to have a model RH problem which resembles the RH problem for P, we let 

It is straightforward to check the RH conditions that are satisfied by 3>. 

RH problem for <3? 

(a) $ : C \ R -> C 2x2 is analytic. 

(b) For i£R, 

o)- (3 - 49) 

(c) <E> behaves as follows as £ — > oo, 

$(() = $(C)C fca3 e- £ ^ <T3 icJi, as ( -> oo, ImA > 0, (3.50) 

= $(C)C fc<J3 e-^ ?£S<T3 , as C ^ oo, ImA < 0, (3.51) 
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where $ has the asymptotic expansion 

8(C) = I + ~£i + ^$2 + 0(C 3 ), (3-52) 

with 

3.5.3 Construction of the parametrix 

We construct P in the form 

P(X) = EkiX^ie-^CiX); k, e- 1/2 s(A))e ± ^( A )-3 _ as ±Im A > 0) (3.55) 
where £ is a real conformal mapping which maps v to 0, s is analytic near v, and 

W = Hf >(^M^-^, «ImA>0. 

l-Pfc ° o) (A)C(A)- fcCT3 e-^ ,J3 , asImA<0, 

which defines analytically near v, with 

E k {v) = {-v) 1/A {v - uy^^e-^NC'ivy^e^M-^iv - u)" fc(T3 . (3.57) 

As before we take 

2 

y = — j — yv — u(x — x), (3.58) 

and let k be the non-negative integer closest to y, and := y — k. Using (|3.4ip . (|3.49|) . 
and (|3.55p . one easily verifies that P satisfies the correct jump condition near v. 

In order to have the right matching (|3.43|) between P and p(°°\ we need to ex- 
ploit the remaining freedom in choosing the conformal mapping C(A) and the analytic 
function s(A) in such a way that 

20(A;x) = -C(A) 2 - 2s(A;x)C(A) - yelne. (3.59) 

If we take a look at 4> in (|3.13|) . we see that 



<^(A; x) = —\/v — u(x — x + ) — (V X — u — \Jv — u){x — x + ) 

+ [\f'L(0 + M)V>^dZ. (3.60) 

J u 

Let us now define £ by 

C(A) 2 = -2^(A;x + ) = -2 [\f L (0 + 6t)y/\^d^, (3.61) 

J u 
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so that 

((v) = 0, C» = (« - u) 1/4 (-#>; u)) 1 / 2 > 0. (3.62) 
Furthermore let 

1 



s(X;x) = (V X — u — y/v — u)(x — x + ), (3.63) 
C[X) 

which defines s analytically near v with 

s{v > x) = 2{ V -um-l{ V] u))^ - (3 - 64) 

Now by ([33DD, HEBIK , and (|5^5]l . it follows that (13391) holds. 

We constructed the parametrix in such a way that for A G 3i7„, 

P(A)P(°°)(A)- 1 = P fe (A) ^ + 1^$! + + 0(e 3 / 2 )^ ^(A)- 1 . (3.65) 

If we choose k to be the non-negative integer closest to y, then it follows from f|3. 56[) 
that P and p(°°) have a good matching at as long as y is not to close to a half 
positive integer. In order to have a uniform matching, also when y approaches a half 
integer, we need to improve the parametrices. 

3.6 Improvement of parametrices 

The goal of this section is to improve the local parametrices and the outside parametrix 
in such a way that they match uniformly for y bounded. This section is inspired by [2]. 
Let us define, if y > —1, 

P(°°)(A) = (l + P (oo) (A), (3.66) 

P U (A) = (l+ C{ ^f ) Pu(\), (3.67) 

where C is a nilpotent matrix (so that the determinant of Q is identically 1) which 
we will determine below. With those improved definitions of the parametrices, p(°°) 
and P u satisfy the same RH conditions as before (see ()3.24|) - (|3.25j) for p(°°)). If C is 
bounded for small e, we have by (|3.38l) 

P u {X)P {oo Xu)- 1 = 7 + 0(eln 2 e)), as A e dU u , (3.68) 

if e — > 0. Next we define the improved local parametrix near v as follows, 

1 — 5e2e s2 



P v (X) = E k (X) | , 22nh2 2 ^ A) 

- ( 1-J )e5e - 2 ^1 1 

x $(e- 1 / 2 C(A);A:,e- 1/2 s(A))e ± ^ (A)<T3 , (3.69) 

as ±ImA > 0, with 

E k (X) = (l + E k (X), (3.70) 
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and with 
S = h 



if k<y<k + ± 5 = 0, \ik-\<y<k. 



(3.71) 



In any case one of the off-diagonal entries of the second factor on the right hand side 
of (|3.69p vanishes. If we let e — > in such a way that y remains bounded, we have that 
s(X;x) = O(elne), and it follows from ([330l - ([33Tjh (pT55]h and (l33H|) that 



P„(A)P (oo) (A)- 1 = / + 0(e 1/2+|Afel ), 



as A G dU v . 



(3.72) 



Note that E k is not analytic at v: it has a simple pole at v. However we can choose the 
matrix C such that P is bounded near v. After a straightforward calculation, it turns 
out that this is the case if 



C = E k {v) 



5d x 
(1 - 5) Cl 



E k (v) 



1 



-(1 - 6)c 2 



-5d 2 
1 



E' k {v) 



with 



ci = ci(A;;eJ ■ 
c 2 = c 2 (k;e) ■ 

d\ = di(k; e) 
d 2 = d 2 (k; e) 



C(v) 
ci("(v) 



2(>(v 
i 



2s(v)s'(v) Cl , 



i s 2 (v) 



2*hlC'(v) 

diC'jv) 
2(>(v) 



+ 2s(v)s'(v)d 1 



5d x 
(1 - 5) Cl 

(3.73) 



(3.74) 
(3.75) 

(3.76) 
(3.77) 



Note that c\{k + 1; e)d\{k\ e) 
particular, writing E = E k , we have 



^77^3 ; and that the matrix C bounded for small e. In 



C 



12 



-(1-5) 



ci-E^u) 



detS(u) - ci [£y i2 (t;)^22(w) - £7 12 (t;)^ 2 (t;)] 



+ 5 



detE(v) - di [£n(tO£ 21 (t>) - ^(^EaiM] 



Observe that 
det#(u) 
E u (vf ■- 

El 2 (vf : 



-2i\/—v, 



v — u 



2k >//■ \— 2fc,— A* 



v — u 



(A(v-u)) 2k (\v) 2k e Ak , 



E u (v)E' 21 (v) - E' u (v)E 21 (v) = 2iV^ - u)) 2fc 1 C»e~ A * 

^ 12 («)^ 2 («) - Sxa^)^^) = 2*^ «(«)(« - n)) 2 "" 1 C'(«)e Afe . 



. (3.78) 

(3.79) 
(3.80) 

(3.81) 

(3.82) 
(3.83) 
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This leads to 



C 12 (k,A k ) = -2{l-S)yfr=H- 17 ^ 



'l + ci7 2fe - 1 C'(f)e Afc 
+ 2SV^ / 1 ] 2k \„,\ a +0(e 2 ln 2 e), (3.84) 

in the double scaling limit where e — > and x — x + = O(elne), with 

7 = 4£ / (u)(v-u). (3.85) 

3.7 Final RH problem 

Now we define R in such a way that it has jumps that are uniformly / + 0(e l l 2 ) in 
the double scaling limit where e — > 0, x — x + = 0(e In e), also when y is close to a half 
integer. Therefore we let 

' S{\-x)P^\\)- 1 , as A G C \ (U u U U v ), 
R(X;x,e) = I S(A; x)P u (A; x)" 1 , as A G U u , (3.86) 
S(X;x)P v (X;x)~ 1 , as A G U v , 

where we use the notations P u = P u and P v = I in the case where y < — 1. Obviously, 
i? is analytic in C \ (Eg U <9C/ U U On <9lT m U the jump matrix for R is close 

to the identity matrix because of the matching of the (improved) local parametrices 
with the outside parametrices. Outside the disks U u and U v , the decay of the jump 
matrix to / is inherited from the decay of the jump matrix for S. Inside the disks, 
there are residual jumps on £5 because the jump matrices for the local parametrices 
are not exactly the same as the ones for S. However, using the small e asymptotics for 
the reflection coefficient, one verifies easily that they are of the form / + 0(e) as e — ► 0. 

RH problem for R 

(a) R is analytic in C \ (£5 U dU u UdU v ). 

(b) R has the jump condition R+(X) = R-(X)vr(X) for A G Eg U dU u U dU v , where 

v R (X) = I + 0(e-i), for \ G Z s \ (U U U U v ), (3.87) 

v R (X) = 1 + 0(e), for A G £ 5 n (U u U C/„), (3.88) 

Wji (A) = / + 0(eln 2 e), for A G dU u , (3.89) 

v R (X) = I + <D(e^ +lAkl ), for AG 9^, (3.90) 

in the double scaling limit where e — ► in such a way that y remains bounded. 
We choose the clockwise orientation for dU u and dU v . 

(c) As A — > 00, we have 

i?(A) = /+^i + 0(A- 2 ). (3.91) 
A 
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All jumps are thus of the form 1 + 0(e 1 ^ 2 ), and if y is a half integer or if y < — 1, 
the situation improves to / + C(eln 2 e). Taking a closer look at the jump for R on dU v , 
and using s(X;x) = O(elne), we obtain using (|3.69p and (|3.52p that 

v R (X) = n(A)p(°°)(A)- 1 



E k {\) \I + 



C(A) \ v <52vr/i|_ 1 

OaFI, 



^TTTwl 2 fc(fc-2) +C(e2lne) E^(A) -1 . 



By (137731) . this leads to 



„ B(A) _ , + g^ (A) r ^ P - ^ + 0(e) . (3.92) 



We can write 

Ufl (A) = / + e5+l A ^t.W(A) + C(eln 2 e), for AG SgU <9C/ U U ffl7 w , 

uniformly for e — > and y bounded. By a standard procedure for small-norm RH 
problems [7], [20], it follows that 12 has a similar expansion in the double scaling limit: 

12(A) = / + e 1 / 2 +l A »=ll2( 1 )(A) + 0(eln 2 e). (3.93) 

Compatibility of (|3,92p with (|3.93p and the jump condition 12+ = R-Vr on dU v gives 
the relation 

R+\X) - R-\ty = Vr\\), fov\edU v . (3.94) 

In addition tfW is analytic in C \ dU v , and R^(oc) = 0. The unique function which 
satisfies those (additive) RH conditions is given by 

R(1) ;v) U\eC\U v , (3.95) 
A — v 

fl(D(A) = g^fe) -vg)(A), S\€U Vi (3.96) 
A — v 

and consequently we have the following asymptotics for 12i defined by (|3.9ip . 

Ri(x,t,e) =eHl A fclRes(vJj 1) ;v) + 0(eln 2 e). (3.97) 

Using the definition (|3.86p of 12, the improved outside parametrix given by (|3,32p . 
and the expansion (]3.66p of the outside parametrix at infinity, we obtain that 

S'ii(A; x, t, e) = 1 — i(Ri,i2(x, t, e) + C12 + 2ky/v — u) = + 0(\~ l ), as A — > 00, 



(3.98) 

which implies by (|3.20p that 

u(x, t,e) = u- 2ed x {Ri,n{x, t, e) + C 12 (x, t, e)) . (3.99) 
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Substituting the asymptotics obtained in (|3.97p gives 

u(x,t,e) = u-2ed x (eS+I^Res^; v) 12 + C\ 2 (x, t, ej) +C(eln 2 e). (3.100) 



It is not so obvious that it is allowed to formally take derivatives of the asymptotics 
obtained for R. However R is analytic as a function of y, and the asymptotic expansion 
(13.931) can be shown to hold true for y in a small complex neighborhood of any i/6t. 
Using this property, (|3.100p can be justified. 
After some calculations we obtain 



2el^.[el Afc lRes(- (1) - 



v R ;v) 12 l 



-8e2 



{v — u) 



1 



2tt^ 7 



- + 2tt^_ i7 



2fc-l 



(3.101) 



where hk is defined in (|3.45p and 7 is defined in (I3,85p , By (13.781) . 



2ed r C 



(v-u)(l-S) 



e/2 



*27T^_ l7 



2fc-l 



1 + e i+A fc27r/l 2_ i7 2fe-l 

8(v - u)S 
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2,,2fc+l> 



(l + ei- A x/(2irhl 1 2k + 1 ] 



(3.102) 



Putting together the above formulas and using the fact that 6 = 1 for k < y < k + ^ 
and (5 = for /c + ^<y</i; + lwe arrive at the expression 

u(x,t,e) = u + 2(v -u) [sech 2 (X fc ) + 4e _2Xft - 1 ] +C(eln 2 e), (3.103) 
\i k < y < k + 1/2, and at the expression 

u(x,t,e) = u + 2(v-u) [sech 2 (X fc ) +4e 2Xfc+1 ] +C(eln 2 e), (3.104) 
if k + 1/2 < y < k + 1, where 

*ifc = \{\ ~ V + *) Ine - ln(v& fc ) - (fc + \) hi 7. (3.105) 

We can write this also in the following more elegant form, which holds for all \y\ < M: 

00 

u(x,t,e) = n + 2(u-n)^sech 2 (X A; ) + C'(eln 2 e). (3.106) 

k=\ 

For each value of y, there are at most two of the sech 2 -terms which are bigger than 
C(eln 2 e), all the others are absorbed by the error term. This proves Theorem 11.21 



Acknowledgements 

TC is a Postdoctoral Fellow of the Fund for Scientific Research - Flanders (Belgium), 
and was also supported by Belgian Interuniversity Attraction Pole P06/02, and by the 
ESF program MISGAM. TG acknowledges support by the ESF program MISGAM. 
TG and TC acknowledge support by ERC Advanced Grant FroMPDEs. 



23 



References 



[1] Beals, R.; Deift, P.; Tomei, C. Direct and inverse scattering on the line, Mathemat- 
ical Surveys and Monographs, 28. American Mathematical Society, Providence, RI, 
1988. xiv+209 pp. ISBN: 0-8218-1530-X 

[2] Bertola, M.; Lee, S.Y. First Colonization of a Spectral Outpost in Random Matrix 
Theory, Constr. Approx. 30 (2009), 225-263. 

[3] Bertola, M.; Tovbis, A. Universality in the profile of the semiclassical limit solu- 
tions to the focusing Nonlinear Schroedinger equation at the first breaking curve, 
larXiv:0909.32641 

[4] Claeys, T. Birth of a cut in unitary random matrix ensembles, Int. Math. Res. 
Not. 2008 (2008), no. 6, Art. ID rnml66. 

[5] Claeys, T.; Grava, T. Universality of the break-up profile for the KdV equation 
in the small dispersion limit using the Riemann-Hilbert approach, Comm. Math. 
Phys. 286 (2009), 979-1009. 

[6] Claeys, T.; Grava, T. Painleve II asymptotics near the leading edge of the oscilla- 
tory zone for the Korteweg-de Vries equation in the small dispersion limit, Comm. 
Pure Appl. Math. 63 (2010), 203-232. 

[7] Deift, P.; Kriecherbauer, T.; McLaughlin, K.T-R ; Venakides, S.; Zhou, X. Strong 
asymptotics of orthogonal polynomials with respect to exponential weights, Comm. 
Pure Appl. Math. 52 (1999), 1491-1552. 

[8] Deift, P.; Venakides, S.; Zhou, X. New result in small dispersion KdV by an 
extension of the steepest descent method for Riemann-Hilbert problems, Int. Math. 
Res. Notices 6 (1997), 285-299. 

[9] Deift, P.; Zhou, X. A steepest descent method for oscillatory Riemann-Hilbert 
problems. Asymptotics for the MKdV equation, Ann. Math. 137 (1993), no. 2, 
295-368. 

[10] Dubrovin, B. On Hamiltonian perturbations of hyperbolic systems of conservation 
laws, II: universality of critical behaviour, Comm. Math. Phys. 267 (2006), no. 1, 
117-139. 

[11] Dubrovin, B. On universality of critical behaviour in Hamiltonian PDEs, Geome- 
try, topology, and mathematical physics 59-109, Amer. Math. Soc. Transl. Ser. 2 
224, Amer. Math. Soc, Providence, RI, 2008. 

[12] Eynard, B. Universal distribution of random matrix eigenvalues near the "birth of 
a cut" transition, J. Stat. Mech. 7 (2006), P07005. 

[13] Faddeev, L.D. Inverse problem of quantum scattering theory II, J. Sov. Math. 5 
(1), 1976, 334-396. 

[14] Fokas, A.S.; Its, A.R.; Kitaev, A.V. The isomonodromy approach to matrix models 
in 2D quantum gravity, Comm. Math. Phys. 147 (1992), 395-430. 



24 



[15] Fujiie, S.; Ramond, T. Matrice de scattering et resonances associees a une orbite 
heterocline (French) [Scattering matrix and resonances associated with a hetero- 
clinic orbit], Ann. Inst. H. Poincare Phys. Theor. 69 (1998), no. 1, 31-82. 

[16] Gardner, S.C.; Greene, J.M.; Kruskal, M.D.; Miura, R.M. Korteweg-de Vries equa- 
tion and generalization, VI. Methods for exact solution, Comm. Pure Appl. Math. 
27 (1974), 97-133. 

[17] Grava, T.; Klein, C. Numerical study of a multiscale expansion of the Korteweg-de 
Vries equation and Painleve II equation, Proc. R. Soc. Lond. Ser. A Math. Phys. 
Eng. Sci. 464 (2008), no. 2091, 733-757. 

[18] Grava, T.; Tian, F.-R. The generation, propagation, and extinction of multiphases 
in the KdV zero-dispersion limit, Comm. Pure Appl. Math. 55 (2002), no. 12, 
1569-1639. 

[19] Gurevich, A. G. ; Pitaevskii, L. P. Non stationary structure of a collisionless shock 
waves, JEPT Letters 17 (1973), 193-195. 

[20] Kamvissis, S.; McLaughlin, K.D.T-R.; Miller, P.D. Semiclassical soliton ensem- 
bles for the focusing nonlinear Schrddinger equation, Ann. Math. Studies 154, 
Princeton Univ. Press, Princeton (2003). 

[21] Lax, P.D.; Levermore, D.D. The small dispersion limit of the Korteweg de Vries 
equation, I,II,III, Comm. Pure Appl. Math. 36 (1983), 253-290, 571-593, 809-830. 

[22] Mo, M.Y. The Riemann-Hilbert approach to double scaling limit of random matrix 
eigenvalues near the "birth of a cut" transition, Int. Math. Res. Not. 2008 (2008), 
no. 13, Art. ID rnn042. 

[23] Ramond, T. Semiclassical study of quantum scattering on the line, Comm. Math. 
Phys. 177 (1996), no. 1, 221-254. 

[24] Tian, F.-R. Oscillations of the zero dispersion limit of the Korteweg-de Vries equa- 
tion, Comm. Pure Appl. Math. 46 (1993), 1093-1129. 

[25] Tsarev, S. P. Poisson brackets and one-dimensional Hamiltonian systems of hy- 
drodynamic type, Dokl. Akad. Nauk. SSSR 282 (1985), 534-537. 

[26] Venakides, S. The Korteweg de Vries equations with small dispersion: higher order 
Lax-Levermore theory, Comm. Pure Appl. Math. 43 (1990), 335-361. 

[27] Whitham, G. B. Linear and nonlinear waves, J. Wiley, New York (1974). 
Tom Claeys 

Cite Scientifique - Laboratoire Painleve M2 
F-59655 Villeneuve d'Ascq, FRANCE 
E-mail : torn . claeysOmath . univ-lillel . f r 

Tamara Grava 
SISSA Via Beirut 2-4 
34014 Trieste, ITALY 
E-mail: grava@sissa.it 



25 



